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A combination of qualitative analysis and numerical study indicates that vacuum 
T 2 symmetric spacetimes are, generically, oscillatory. 

Thirty years ago, Bclinskii, Khalatnikov and Lifshitz proposed that the dynam- 
ics of spatially inhomogeneous solutions to Einstein's equation near a cosmological 
singularity are oscillatory in general. H The quantities which oscillate arc the gen- 
eralized Kasner exponents, which are the eigenvalues of the extrinsic curvature 
divided by the mean curvature. In various known cases, there exists a foliation and 
threading of a solution to Einstein's equation in the neighborhood of a cosmologi- 
cal singularity such that the generalized Kasner exponents along each thread either 
converge to a limit at the singularity (if this occurs on all threads we call the solu- 
tion convergent), or approximately follow the BKL sequenced Along threads with 
oscillatory behavior, the BKL sequence is realized because the evolution becomes 
approximately Kasner, but then there is inevitably a transition to a different Kasner 
evolution. Each transition can be approximated, yielding the BKL sequence. The 
singularity is at finite proper time, but there are an infinite number of oscillations. 

Numerical simulations of vacuum T 2 symmetric spacetimes show that an 
"asymptotic regime" is reached in which the generalized Kasner exponents fol- 
low a portion of the BKL sequence. Qualitative analysis in the asymptotic regime 
indicates that the oscillations will, in general, continue without, &ud at almost every 
spatial point. These methods have been described previously.BBa 

Gowdy spacetimes are vacuum T 2 symmetric spacetimes such that the (space- 
like) symmetry orbits are surface orthogonal. It is thought that all Gowdy space- 
times are convergent. We find oscillatory dynamics in non-polarized T 2 symmetric 
spacetimes if the symmetry orbits are not surface orthogonal, in which case the 
spatial topology musLbe the three torus. A global foliation and threading of these 
spacetimes is known.cl The metric can be written 

g = -e^dr 2 + e h±J ^de 2 +ae p - T [dx + Qdy+{G 1 + QG 2 )dO 

- {Ah + Q M 2 ) e~ T dr] 2 + a e- p - T [dy + G 2 dO - M 2 e~ T dr] 2 , ' ' ' 

with the singularity in the direction of increasing r. The metric functions are 
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independent of x and y and periodic in (L If Q vanishes everywhere the solution is 
polarized, and thought to be convergent.!] If the non-positive function d T fi vanishes 
everywhere, the symmetry orbits are surface orthogonal. If \d T /j.\ <C 1, then the 
generalized Kasner exponents are approximately 

d T X + d T fi+l 2{v-l) -2(v + l) 

d T \ + d T fi — 3 drX + dr^ — i d T \ + d T (i — 3' 

The error is less than 30^/— 9 r /i, which we show studying perturbations of linear 
operators. The denominator in (||) is strictly negative. 

Numerical simulations of a spacetime with metric (Q) show that the evolution 
becomes approximately Kasner at each value of 8. The signature for this is that 
both /i and also v = y/ '{d T P) 2 + e 2P (d T Q) 2 are approximately constant in time, and 
d T \ rs — v 2 . Both numerical simulations and qualitative analysis indicate that there 
are two types of transition. The signature of a transition is that v is not constant 
in time and d T X is not approximated by — v 2 . For each type of transition a rule for 
v is obtained by considering an exact "transition solution" which approximates the 
evolution. One type of transition is driven by the spatial curvature. It has been 
studied in various classes of spacetimes, including the Gowdy spacetimes, and gives 
two consecutive Kasner spacetimes in the BKL sequence. In the Gowdy spacetimes 
the oscillations eventually cease. If d T n ^ then a second type of transition 
occurs, in which d T \x grows in magnitude and decays again, and throughout which 
the evolution remains Kasner, but the Kasner directions rotate. This rotation is 
geometrical. A given Kasner direction may be tangent to the symmetry orbit before 
the transition, but not tangent after. In the spatially .homogeneous setting this type 
of transition was noticed using Hamiltonian methoddj ( "centrifugal potential" ) and 
has since been noticed using dynamical systems methods.tl 
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